Stability of the Fast Multipole Method for Helmholtz Equation in Three Dimensions
نویسنده
چکیده
Stability limits for the diagonal forms approximating the free space Green’s function in Helmholtz’ equation are derived. It is shown that while the original approximation of the Green’s function is stable except for overflows, the diagonalized form becomes unstable due to errors from roundoff, interpolation, choice of quadrature rule and approximation of the translation operator. Numerical experiments validate the theoretical estimates.
منابع مشابه
Remarks on the implementation of the wideband FMM for the Helmholtz equation in two dimensions
This paper describes a simple version of the Fast Multipole Method (FMM) for the Helmholtz equation in two dimensions. We discuss both the underlying theory and some of the practical aspects of its implementation to allow for stability and high accuracy at all wavelengths.
متن کاملA wideband fast multipole method for the Helmholtz equation in three dimensions
We describe a wideband version of the Fast Multipole Method for the Helmholtz equation in three dimensions. It unifies previously existing versions of the FMM for high and low frequencies into an algorithm which is accurate and efficient for any frequency, having a CPU time of O(N) if low-frequency computations dominate, or O(N logN) if high-frequency computations dominate. The performance of t...
متن کاملAn Implementation of Low-Frequency Fast Multipole BIEM for Helmholtz’ Equation on GPU
Acceleration of the fast multipole method (FMM), which is the fast and approximate algorithm to compute the pairwise interactions among many bodies, with graphics processing units (GPUs) has been investigated for the last couple of years. In view of the type of kernel functions, the non-oscillatory kernels (especially, the Laplace kernel) were studied by many researchers (e.g. Gumerov), and the...
متن کاملCubic spline Numerov type approach for solution of Helmholtz equation
We have developed a three level implicit method for solution of the Helmholtz equation. Using the cubic spline in space and finite difference in time directions. The approach has been modied to drive Numerov type nite difference method. The method yield the tri-diagonal linear system of algebraic equations which can be solved by using a tri-diagonal solver. Stability and error estimation of the...
متن کاملA broadband fast multipole accelerated boundary element method for the three dimensional Helmholtz equation.
The development of a fast multipole method (FMM) accelerated iterative solution of the boundary element method (BEM) for the Helmholtz equations in three dimensions is described. The FMM for the Helmholtz equation is significantly different for problems with low and high kD (where k is the wavenumber and D the domain size), and for large problems the method must be switched between levels of th...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
عنوان ژورنال:
دوره شماره
صفحات -
تاریخ انتشار 2003